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Abstract. Speckle interferometric observations made with the 6 m telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciences in 2000 revealed the triple nature of the nearby {iiHip = 
51.80 ± 1.74 mas) low-mass young (« 200 Myr) star GJ900. The configuration of the triple system al- 
lowed it to be dynamically unstable. Differential photometry performed from 2000 through 2004 yielded I- 
and K-haxiA absolute magnitudes and spectral types for the components to be /a=6.66±0.08, /b=9.15±0.11, 
/c=10.08±0.26, ii'A=4.84±0.08, A"s=6.76±0.20, A'c=7.39±0.31, S'pa~K5-K7, SpB^MZ-MA, Spc~M^-MQ. 
The "mass-luminosity" relation is used to estimate the individual masses of the components: A1a~O.64A10, 
Mb~Q.'2^Mq, A1c~O.13Al0. From the observations of the components relative motion in the period 2000- 
2006, we conclude that GJ 900 is a hierarchical triple star with the possible orbital periods Pa-sc~80 yrs and 
P_Bc~20 yrs. An analysis of the 2MASS images of the region around GJ 900 leads us to suggest that the system 
can include other very-low-mass components. 



1. INTRODUCTION 

In the solar neighborhood M-type dwarfs make up about 
70% and no less than 40% of all population in terms of 
the number of stars and mass, respectively. The lifetime of 
the lowest-luminosity main-sequence stars exceeds the age 
of the Universe, making them good candidate objects for 
the study of the properties of the Galactic disk including 
the history of star formation in the local volume. 

The study of young M-type field dwarfs may lead 
to the discovery of systems with brown dwarfs. The in- 
terest toward the discovery and study of multiple sys- 
tems with substellar components has been growing pro- 
gressively since the discovery of the first brown dwarf 
GJ229 B dNakajima et al., 1995| . It is evident that such 
objects can be found only at small distances from the 
Sun. An example of such objects is the group of brown 
dwarfs GJ569 B, which had its first orbit and dynamical 
masses estimated from the results of observations made 
with three major telescopes: the 6 m Bolshoi Azimuthal 
Telescope (BTA) of the Special Astrophysical Observatory 
of the Russian Academy of Sciences, Keck II, and MMT 
dKenworthy et al., 200T| ). 

Until recently, the review by Henry and McCarthy 
( [Henry fc McCarthy, 1993| remained the most repre- 
sentative summary of empirical data on the masses 
and luminosities of cool dwarfs obtained using various 
observational techniques. In recent years, new data on 



the principal parameters of low-mass stars in binary 
and multiple systems have been published. These data 
were obtained by combining different observational meth- 
ods: adaptive optics imaging combined with accurate 
radial- velocity measurements ( [Segransan et al., 2000 ) , 
speckle and long-baseline interferometry with ra- 
dial velocities measurements ( |Balega et al., 2007| 
[Boden et al., 2006| , and space astrometry performed 
with fine guidance sensors of the Hubble Space Telescope 
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( [Torres et al, 1999| [Bmedict et al., 2000| ). These studies 
allowed the main empirical relations to be substan- 
tially refined Pelfosse et al., 2000D . This concerns, m 
particular, the "mass-luminosity" relation, which is of 
great importance for the investigation of parameters 
of individual stars and of the entire population of our 
Galaxy. 

In 1998, a speckle interferometric survey of low-mass 
binaries and suspected binaries discovered by Hipparcos 
( [Ferryman et al., 1997[ ) astrometric satellite was started 
at the BTA 6 m telescope ( [Balega et al., 2002[ ). 
Observations were performed at visual and infrared wave- 
lengths with the aim to identify pairs with fast relative 
motion of components, which allow model-independent 
masses to be determined over a short time interval. In 
addition to measuring relative component positions with 
an accuracy of about 1-2 milliarcseconds (mas), we also 
measured the magnitude differences in the V, R, I, J, 
H, and K bands for most of the binaries. Objects of 
this program include the red star G J 900=Hip 116384 
with the Hipparcos ( [Ferryman et al., 1997[ ) parallax of 



2 



Malogolovets et al.: GJ900: A new hierarchical system with low-mass components 



7rmp=51.80±1.74 mas. The results of astrometric mea- 
surements of this star allowed it to be suspected as a bi- 
nary based on a number of indicators (the object has flag 
'S' in the Hipparcos catalog; see ( [Lindegren, 1997P ) and 
that is why we included it into our program list. 

The very first interferometric observations performed 
in November 2000 with BTA showed GJ 900 to be a triple 
system ( Balega et al., 2006 ). On the reconstructed image 
two fainter components are located 0.5" and 0.7" from 
the main star. The compact configuration of the GJ 900 
system may indicate that it belongs to the class of dynam- 
ically unstable multiple stars. One must, however, bear in 
mind that the components of this system may appear to be 
located at a comparable distances from the main star only 
in the sky-plane projection. To verify these assumptions, 
we included this system into the program of monitoring 
of the relative motion of components. 

In 2002-2003, Martin ( [Martin, 200^ observed GJ 900 
with the CIAO adaptive optics system of the 8.2 m Subaru 
Telescope and confirmed the two faint companions of the 
central object. The results of two observations in the in- 
frared bands H and K separated only by a five-month 
interval led him to conclude that the system is dynami- 
cally bound. 

In this paper we report the results of the interferomet- 
ric measurements of positional parameters and magnitude 
differences of the components of G J 900 made during the 
period from November, 2000 through December, 2006, de- 
termine the absolute magnitudes and estimate the masses 
of the stars. We analyze the possible dynamical stability 
of the system based on the observations performed. 

2. OBSERVATIONS AND DATA ANALYSIS 

Speckle interferometric observations of GJ 900 were made 
with the BTA 6 m telescope in the V (A/AA=550/30 
nm), / (A/AA=800/100 nm), and A'(A/AA=2115/214 
nm) bands, where A is the central wavelenght, AA is the 
half- width of the band. A fast 512x512 Sony ICX085 
CCD combined with a three-stage image intensifier was 
used as a detector from 2000 through 2004. During our 
2006 observations we employed a new EMCCD (Electron 
Multiplying Charge Coupled Device) system with higher 
quantum efficiency and linearity. The image scale was 
equal to 4.1 and 6.7 mas/pixel for the first and second fa- 
cilities, respectively. The exposure time of speckle interfer- 
ograms varied from 5 to 20 msec depending on brightness 
of an object and seeing conditions. Infrared observations 
were made with HAWAII infrared detector of the Max- 
Planck- Institute for Radio Astronomy. Table[T] gives a log 
of observations. This table gives the following data for each 
measurement: date as a fraction of Besselian year, seeing 
P in arcseconds, the number of speckle interferograms in 
each series, and the name of the filter. 

The angular distances p, position angles 0, and magni- 
tude differences between the components Am from speckle 
interferometric observation with the BTA telescope are 
given in Table 2. Due to low signal-to-noise ratio in the 
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Fig. 1. Images of GJ900 reconstructed from the 
2115/214-nm filter observations with the BTA 6 m tele- 
scope in August 2003 (left) and from the 800/100-nm in 
December 2006 (right). North is at the top and East on 
the left. 



2006.9465 measurements in T^-band, the relative positions 
for this data are not presented in the table. A descrip- 
tion of the technique of the determination of relative posi- 
tions and component magnitude differences inferred from 
the averaged over the series power spectra of speckle in- 
terferograms can be found in the paper by Balega et al. 
QBalega et al., 2002] ). The resolution diffraction limit was 
equal to 0.022", 0.033", and 0.088" in the V, I, and K 
bands, respectively. The accuracy of the measurement of 
position parameters is equal to 0.3-1.0° and 3-8 mas in 
position angle 9 and angular separation p, respectively. 
The errors of measured 9 and p depend on a number 
of parameters: component separation, magnitude differ- 
ences, and seeing p. The accuracy of the determination of 
component magnitude differences from the reconstructed 
power spectra is also a function of the same parame- 
ters. This accuracy varies from 0.05 to 0.2 for objects 
with mv=S-10. The modulus of the Fourier transform 
of the object (visibility) was obtained from the series of 
speckle interferograms with the classical speckle intcrfer- 
ometry method. For image reconstruction we used the 
bispectrum speckle interferometry method ( [Weigelt, 1977[ 
[Lohmann et al., 1983D . Fig. [T] shows the reconstructed in- 
frared and visual images of GJ 900 based on observations 
made in 2003 and 2006. 

To specify the spectral type of the object, the 3600 - 
6200 A spectrum of the object was taken in October 2006 
with the UAGS spectrograph of the Zeiss-1000 telescope 
with a dispersion of 1.35 A/pixel [2l 
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Fig. 2. Spectrum of GJ 900 in the 3600 - 6200 A wavelength interval taken with the UAGS spectrograph of Zeiss-1000 
telescope in October 2006. The strongest absorption Hnes and bands — those of TiO, MgH, Cal, and Nal — are 
indicated. 



3. ABSOLUTE MAGNITUDES, MASSES AND 
SPECTRAL TYPES OF THE COMPONENTS 

To compute the /- and X-band absolute magnitudes of 
the components of GJ 900, we used the results of differ- 
ential speckle photometry performed with the BTA tele- 
scope (Table 2) and published integrated magnitudes of 
the system in these bands. For the sake of complete- 
ness, this table also includes the photometric data from 
dMartin, 20"03| . The component magnitude differences de- 
termined using the results of speckle interferometry per- 
formed with BTA agree, on the whole, with those inferred 
from observations performed with the Subaru Telescope 
adaptive optics. At the same time, the results of differ- 
ential photometry do not rule out the variability of one 
or several components of the system. The integrated I- 
band magnitude determined from the color index of the 
system, V — I = 1.65 ([Bessel, 1990[), and the visual mag- 



nitude adopted from ( Salim fc Gould, 2003 | is equal to 
7717=7.94. The integrated if-band magnitude is equal to 

=6. 01 ±0.01 dAlonso et al., 1994D . We combine these 
data with the Hipparcos parallax and magnitude differ- 
ences from Table 2 to infer the following estimates for the 
absolute magnitudes of the components of GJ 900: 

/a=6.66±0.08, 7^a=4.84±0.08, 

/i3=9.15±0.11, Ki3=6.76±0.20, 

/c=10.08±0.26, Kc=7.39±0.31, 

where for the / band we used the mean result aver- 
aged over two observations made with the BTA 6 m tele- 
scope: Amf ^=2.49±0.07, and Amf^ = 3.42 ± 0.24. The 



i?-band component magnitude differences measured using 
adaptive optics ( [Martin, 2003D yield the following absolute 
magnitudes for the components: 

iJA=5.11±0.08, 

iJs=6.85±0.09, 

iJc=7.54±0.15. 

The luminosities of the components allow their masses 
to be estimated from the "mass-luminosity" relation, how- 
ever, to do this, we must know the age and metallicity of 
the system. 



Gizis et al. (Gizis et al., 2002) compared the activ- 



ity of a large sample of nearby M-type field dwarfs with 
the activity of M-type dwarfs in open clusters and cali- 
brated the "age-activity" relation. We use the V ~ I — 
—6.91 -I- 1.05log{T) relation to infer an age of lO'* Myr for 
the adopted distance modulus of m-A/=1.43. X-ray lumi- 
nosity of late-type stars is yet another indicator of their 
activity. According to ROSAT observations, the X-ray lu- 
minosity of GJ900 is equal to = 108.8 x 10^'' erg/s 
( [Huensch et al., 1999p . This value corresponds to the M- 
dwarf age of 100-200 Myr. 

The fluxes in the Call H and K and Mgll h and k 
lines are good indicators of stellar activity and hence 
of stellar age [Soderblom et al., 1991| [Wilson, 19631 ). 
However, no relations for age as a function of these 
indicators have been derived for stars of spectral types 
later than K. The only solution is to compare GJ 900 
with stars having similar spectra with the age deter- 
mined using other methods. For GJ900 Giampapa et 
al. ( [Giampapa et al., 1989[ ) report Call H and K line 
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Table 2. Speckle interferometric and adaptive optics measurements of GJ900 
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vector 


ir^s 


m^s 


61° 




Am 
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Reference 


2000 8754 


AB 
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3 


316.4 


0.4 


2.42 


0.15 


/ 


(|Balega et al., 2006 p 




AC 
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5 


344.3 


0.4 


3.65 


0.22 






BC 
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6 


13.6 


0.6 










2002 5QQn 


AB 


510 


10 


324.5 


0.1 


1.78 


0.02 




(IMartin 200311 




AC 


760 


10 


344.0 


0.1 


2.55 


0.03 
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AB 
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10 


324.5 


0.1 


1.61 


0.03 


K 


(IMartin 20031) 




AC 
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0.04 
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AB 
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20 


327.4 


0.1 


1.70 


0.04 




(IMartin 200311 
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20 
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0.06 
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AB 
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5 
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1.92 


0.18 
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AC 
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BC 
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4 
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1.0 


0.63 


0.35 
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AB 
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4 


331.7 


0.4 






J 






AC 
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5 


345.1 


0.4 












BC 


224 


6 


20.4 


0.6 










2004.8208 


AB 


606 


3 


335.8 


0.3 


2.56 


0.06 


I 


This paper 




AC 


714 


4 


345.5 


0.4 


3.18 


0.22 








BC 


155 


5 


26.7 


0.5 










2006.9465 


AB 


751 


3 


342.5 


0.3 






I 


This paper 




AC 


708 


8 


344.7 


0.7 












BC 


51 


9 


130.3 


0.8 











flux measurements: log{Fcaii) = 5.18 erg/sxcm^. We 
adopted the Mgll-line flux log{FMgii) = 5.83 erg/sxcm^ 
from dPanagi &: Mathioudakis, 1993D . For compari- 
son, we selected two stars of similar spectral types: 
GJ212 and GJ879. The Call H and K line fluxes of 
these stars are equal to log{Fcaii) = 5.49 erg/sxcm^ 
( [Stauffer fc Hartmann, 1986} [Rutten et al., 1989[) and 
log{Fcaii) = 5.97erg/sxcn? ( [Robinson et al., 1990[ 
[Rutten et al., 1991D , respectively. The Mgll h and 
k line fluxes are equal to log{FMgii) = 5.36 
erg/sxcm^ ( Panagi fc Mathioudakis, 1993[ ) and 
log{FMgii) = 5.98erg/sxcm^, respectively 

( [Panagi k, Mathioudakis, 1993[ ). The ages of these 
stars were inferred from kinematics, isochrones, and 
lithium abundance and vary from 100 to 200 Myr. 



Zuckerman et al. ( [Zuckerman et al., 2006P believe 
GJ 900 to be a possible member of the Carina-Near mov- 
ing group with an estimated age of 200 ± 50 Myr. 
However, the star's membership in this group is doubtable. 
The radial velocity of GJ900 is equal to -10 km/s and 
differs strongly from that of the core of the group (-1-20 
km/s ( [Zuckerman et al., 2006[ )). The center of the moving 
group is at a distance of 30-50 pc from us, whereas GJ 900 
is located within mere 19 pc from the Sun. The equivalent 



width of the Lil A6708 line ( [Zboril fc Byrne, 1998[ ) in the 
spectrum of GJ 900 is several factors of ten smaller than in 
that of the corresponding lines in the spectra of the main 
members of the group. 



Martin ( [Martin, 2003P used the results of the kinemat- 
ical survey of Montes et al. ( [Montes et al., 200 1[ ) to esti- 
mate the age of the star at 50-100 Myr. 



Thus all the available observational data indicate that 
GJ900 is a young system with the age of about 200±100 
Myr. 

Since G J 900 is located in the immediate proximity 
to the Sun (d=19.3 pc) and belongs to the galactic disk 
population, we can suppose that its metal abundance is 
close to the Sun's value. This assumption is supported by 
the results of spectroscopic study of the system performed 
by Zboril and Byrne (|Zboril fc Byrne, 1998). They give 



the metallicity of the star [M/H]=-0.1±0.2. 



Zboril and Byrne (Zboril & Byrne, 1998) determined 



the effective temperature of GJ900 from sensitive pho- 
tospheric lines and molecular bands with the allowance 
for surface gravity and microturbulence. Their esti- 
mate, Te//=4000K, is lower by 200 K than the tem- 
perature inferred from the B-V and R-I color indices 
( [Gliese fc Jahreiss, 1991[ ). It is evident that this tempera- 
ture and the corresponding spectral type K7 refer to the 
main component of the system. 

Speckle interferometric measurements of the /— , K- 
band component magnitude differences and adaptive op- 
tics measurements in the H and K bands allow the ef- 
fective temperature and hence the spectral type of the 
primary component to be estimated using the calibrated 
relation between temperature and V-I and V-K color in- 
dices ( [Alonso et al., 19^6] ). To estimate the T^-band lumi- 
nosity of the primary component, we used the /, H, and 
if-band absolute magnitudes and theoretical isochrones 
for the age of 200 Myr ( [Baraffe et al., 1998[ ). The mean 
Va inferred from three isochrones is equal to 8.35±0.07. 
In this case, given the distance modulus of m — M = 1.43, 
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Fig. 3. Measurement of the angular separation and position angle of GJ 900 B relative of GJ 900 A. Circles and 
asterisks show the speckle interferometric measurements made with the BTA 6 m telescope and the measurements 
made with the Subaru Telescope ( [Martin, 2003^ using adaptive optics, respectively. The bars show the measurement 
errors. 
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Fig. 4. Same as Fig. 3, but for GJ 900 C. 



the color indices of the primary should be equal to {V~ 
J)^=1.70±0.08 and {V-K)A=S.52±0m. The calibration 
of temperature in terms of color index derived by Alonso 



et al. ( [Alonso et al., 1996[ ) depends only slightly on the 
metallicity of the star. We assume that iron abundance is 
equal to the solar value to infer from the estimated {V- 
I)a and {V-K)a indices the temperature of component 
GJ 900 A, T^^.^=4079±180 K, which corresponds to a late 
K-type dwarf. 

An analysis of the GJ 900 spectrum obtained with the 
UAGS spectrograph on the Zeiss-1000 telescope showed 
that the energy distribution and relative intensities of sin- 
gle strong Hnes correspond to those of a K5-K7-type star 
(Fig. [2|). The characteristic features of the spectrum in- 
clude strong TiO and MgH absorption bands. 

For the age of 100-200 Myr and solar chemical 
composition the Mi — and Mk — M. evolutionary 
tracks for low-mass stars computed by Baraffe et al. 
( [Baraffe et al., 1998[ ) imply a primary-component mass in 
the interval from 0.64 to 0.67 Mq. This mass agrees best 



with the primary spectral type of K5 - K7. The masses of 
lower-mass components B and C inferred from the same 
tracks show a much greater scatter. The mass of GJ900 
B is estimated to range from 0.28 to 0.34 M.q, and that 
of GJ900 C, from 0.16 to 0.24 Mq. Note the scatter of 
possible masses is even greater when estimated from the 
K-hasvA photometry of the stars studied. The resulting 
mass estimates correspond to the spectral types of M3 - 
M4 and M5 - M6 for the second and third components, 
respectively. 



4. RELATIVE MOTION OF COMPONENTS, 
LIKELY ORBITAL PERIODS, AND THE 
DYNAMICAL STABILITY OF THE SYSTEM 

The proper motion of GJ900 is equal to 344 mas/yr. If 
we were dealing with accidental projection, components 
B and C would have shifted by almost 2" relative to 
component A during the observing period from 2000 to 
2006. However, the mutual positions of these components 
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Fig. 5. J, H, and if -band 2MASS images of the 2.1'x2.1' region around GJ900 ( |Skrutskie et al., 2006[ ). North is at 
the top and East on the left. 



changed insignificantly during our monitoring program 
(Fig. [3] and 2]). The average annual variation in the posi- 
tional parameters of component B relative to component 
A was equal to 4.3° and 55 mas in position angle and 
angular separation, respectively. The variation of the po- 
sition of component C relative to component B is equal 
to 19.2°/yr and 57 mas/yr for 6 and p, respectively. It 
follows from the relative position changes that the com- 
ponents B and C form the inner short-period subsystem, 
which moves with the component A around the GJ 900 
mass center. The average annual variations of the posi- 
tional parameters imply an orbital period of about 80 and 
20 years for the subsystems A-BC and BC, respectively. 
Thus GJ 900 is a gravitationally bound hierarchical mul- 
tiple system. The planes of the orbital motion of compo- 
nents in the subsystem BC and of component A are tilted 
most probably at a large angle to each other, resulting in 
the observed configuration. 

5. THE PRESENCE OF EXTRA COMPONENTS 
IN THE SYSTEM 



To find the eventual faint 
GJ900 system, we analyzed 
dSkrutskie et al., 2006[) taken 



components in the 
the 2MASS images 

in August, 2000. We 

found on the J-, and isT-band images a faint compan- 
ion of Ril2-13 magnitude ~ 12" Northeast of the central 
object (Figl5]). The K-h&nd image also shows another 
component, at « 15" South of the central source. The 
probability for a star to be located accidentally within 
the 30-arcsec field in the region studied is equal to about 
one percent. Hence GJ 900 is very likely to be a quadruple 
or even a quintuple system. From the intensity ratio it 
follows that the faint components are late M dwarfs. 

To verify these hypotheses, /-band images of the 
GJ 900 vicinity were taken by A.V. Moiseev in February 
2007 with the SCORPIO focal reducer on the BTA tele- 
scope. We took and then averaged 25 ten-second expo- 
sures. The limiting magnitude of the resulting image is 
17™, but we found no objects at the location of the 
Northeastern component. This result suggests that either 



this component is too red to be seen in the / band — 
and this assumption is supported by the absence of the 
object on the /-band plate of the DSS2 survey — or that, 
because of the proper motion of the component or that of 
GJ 900, this object is now projected onto GJ 900. Hence 
infrared photometry is needed to make the final conclusion 
about the membership of components in the multiple sys- 
tem. If the components found do not belong to the multi- 
ple system G J 900, they should shift considerably relative 
to GJ 900 since the time when the corresponding 2MASS 
images were taken. If the relative positions of these com- 
ponents remain unchanged, this object would be a unique 
low-mass multiple system, which is of interest for testing 
the theory of star formation and dynamical evolution of 
stars. 



6. CONCLUSIONS 

Speckle interferometric observations made with BTA 
6 m telescope during the period from 2000 to 
2006 showed that GJ 900 is a gravitationally bound 
triple star. This system belongs to the population 
of the thin disk of the Galaxy and has an age of 
200±100 Myr. The absolute magnitudes of the compo- 
nents are: /a=6.66±0.08, /b=9.15±0.11, /c=10.08±0.26, 
Ka=4.M±Q.08, Kb^6.76±0.20, and /Cc=7.39±0.31 and 
correspond to the spectral types of S'pA~K5-K7, 
SpB-MS-M-i, and Spc~Mb-M 6. We used the ev olu- 
tionary tracks of Baraffe et al. (Baraffe et al., 1998) to 
compute the component masses for the solar metallic- 
ity and for the age of 200±100 Myr: Ma'=^OM-0. 67 Mq, 
Mb~0.28-0MMq, Mc~0.16-0.24Mq. The estimated 



masses and absolute magnitudes of the components agree 
with the results obtained in the H and K bands with the 
Subaru Telescope adaptive optics ( [Martin, 2003p . 

We estimated the orbital periods of the components. 
Pec ~20 yr and Pa-bc ~80 yr, and concluded that 
GJ 900 is an hierarchical multiple star. We explain the 
comparable mutual angular separations between the com- 
ponents by the effect of sky-plane projection. 



Malogolovets et al.: GJ900: A new hierarchical system with low-mass components 



7 



We examined the J-, H-, and -fC-band images of the 
2MASS survey and found possible distant components of 
the GJ 900 system at 12" and 15". If these components are 
gravitationally bound to the triple star, GJ 900 should be 
a young quintuple system of red dwarfs. 
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